Hückel Molecular Orbital Analysis for Stability and Instability of Stacked Aromatic and Stacked Antiaromatic Systems.

Face-to-face stacking of aromatic compounds leads to stacked antiaromaticity, while that of antiaromatic compounds leads to stacked aromaticity. This is a prediction with a long history; in the late 2000s, the prediction was confirmed by high-precision quantum chemical calculations, and finally, in 2016, a π-conjugated system with stacked aromaticity was synthesized. Several variations have since been reported, but essentially, they are all the same molecule. To realize stacked aromaticity in a completely new and different molecular system and to trigger an extension of the concept of stacked aromaticity, it is important to understand the origin of stacked aromaticity. The Hückel method, which has been successful in giving qualitatively correct results for π-conjugated systems despite its bold assumptions, is well suited for the analysis of stacked aromaticity. We use this method to model the face-to-face stacking systems of benzene and cyclobutadiene molecules and discuss their stacked antiaromaticity and stacked aromaticity on the basis of their π-electron energies. By further developing the discussion, we search for clues to realize stacked aromaticity in synthesizable molecular systems.

Exploring Metal Nanocluster Catalysts for Ammonia Synthesis Using Informatics Methods: A Concerted Effort of Bayesian Optimization, Swarm Intelligence, and First-Principles Computation.

This paper details the use of computational and informatics methods to design metal nanocluster catalysts for efficient ammonia synthesis. Three main problems are tackled: defining a measure of catalytic activity, choosing the best candidate from a large number of possibilities, and identifying the thermodynamically stable cluster catalyst structure. First-principles calculations, Bayesian optimization, and particle swarm optimization are used to obtain a Ti8 nanocluster as a catalyst candidate. The N2 adsorption structure on Ti8 indicates substantial activation of the N2 molecule, while the NH3 adsorption structure suggests that NH3 is likely to undergo easy desorption. The study also reveals several cluster catalyst candidates that break the general trade-off that surfaces that strongly adsorb reactants also strongly adsorb products.

Complexation-Triggered Fluctuation of π-Conjugation on an Antiaromatic Dicyanoanthracene Dianion.

The formation of Lewis pairs is an important chemical concept. Recently, the complexation of Lewis acidic tris(pentafluorophenyl)borane with Lewis basic moieties and subsequent reduction has emerged as a fascinating strategy for designing novel reactions and structures. The impact of the complexation and subsequent reduction of antiaromatic systems bearing Lewis base moieties has been investigated. We found how Lewis adduct formation stabilizes an antiaromatic system consisting of 9,10-dicyanoanthracene and tris(pentafluorophenyl)borane by using synthesis, X-ray crystallography, spectroscopic analysis, and quantum chemical calculations.

Frontier Orbital Views of Stacked Aromaticity.

Recent studies have theoretically and experimentally demonstrated that antiaromatic molecules with 4n π electrons exhibit stacked aromaticity according to π-π stacking when arranged in a face-to-face manner. However, the mechanism of its occurrence has not been clearly studied. In this study, we investigated the mechanism of stacked aromaticity using cyclobutadiene. When the antiaromatic molecules are stacked in a face-to-face manner, the orbital interactions between the degenerate singly occupied molecular orbitals (SOMOs) of the monomer unit cause a larger energy gap between the degenerate highest-occupied molecular orbitals (HOMOs) and the lowest-unoccupied molecular orbitals (LUMOs) of the dimer. However, the antiaromatic molecules are more stable in less symmetric conformations, mainly because of pseudo-Jahn-Teller distortions. In the case of cyclobutadiene, the two SOMOs of the monomer unit split into HOMO and LUMO because of the bond alternation. When the molecules are stacked in a face-to-face manner, the HOMO-LUMO gap of the dimer is smaller than that of the monomer due to the interactions between the HOMOs and LUMOs of the two monomer units. When the monomer units are within a specific distance of each other, the HOMO and LUMO of the dimer, which correspond to antibonding and bonding between the units, respectively, are interchanged. This alternation of molecular orbitals may result in an increase in the bond strength between the monomer units, exhibiting stacked aromaticity. We demonstrated that it is possible to control the distance exhibited by stacked aromaticity by engineering the HOMO-LUMO gap of the monomer units.

Molecular Understanding of Adhesion of Epoxy Resin to Graphene and Graphene Oxide Surfaces in Terms of Orbital Interactions.

The adhesion mechanism of epoxy resin (ER) cured material consisting of diglycidyl ether of bisphenol A (DGEBA) and 4,4'-diaminodiphenyl sulfone (DDS) to pristine graphene and graphene oxide (GO) surfaces is investigated on the basis of first-principles density functional theory (DFT) with dispersion correction. Graphene is often used as a reinforcing filler incorporated into ER polymer matrices. The adhesion strength is significantly improved by using GO obtained by the oxidation of graphene. The interfacial interactions at the ER/graphene and ER/GO interfaces were analyzed to clarify the origin of this adhesion. The contribution of dispersion interaction to the adhesive stress at the two interfaces is almost identical. In contrast, the DFT energy contribution is found to be more significant at the ER/GO interface. Crystal orbital Hamiltonian population (COHP) analysis suggests the existence of hydrogen bonding (H-bonding) between the hydroxyl, epoxide, amine, and sulfonyl groups of the ER cured with DDS and the hydroxyl groups of the GO surface, in addition to the OH-π interaction between the benzene rings of ER and the hydroxyl groups of the GO surface. The H-bond has a large orbital interaction energy, which is found to contribute significantly to the adhesive strength at the ER/GO interface. The overall interaction at the ER/graphene is much weaker due to antibonding type interactions just below the Fermi level. This finding indicates that only dispersion interaction is significant when ER is adsorbed on the graphene surface.

Extraordinary Acceleration of an Electrophilic Reaction Driven by the Polar Surface of 2D Aluminosilicate Nanosheets.

To increase chemical reaction rates, general solutions include increasing the concentration/temperature and introducing catalysts. In this study, the rate constant of an electrophilic metal coordination reaction is accelerated 23-fold on the surface of layered aluminosilicate (LAS), where the reaction substrate (ligand molecule) induces dielectric polarization owing to the polar and anionic surface. According to the Arrhenius plot, the frequency factor (A) is increased by almost three orders of magnitude on the surface. This leads to the conclusion that the collision efficiency between the ligands and metal ions is enhanced on the surface due to the dielectric polarization. This is surprising because one side of the ligand is obscured by the surface, so the collision efficiency is expected to be decreased. This unique method to accelerate the chemical reaction is expected to expand the range of utilization of LASs, which are chemically inert, abundant, and environmentally friendly. The concept is also applicable to other metal oxides which have polar surfaces, which will be useful for various chemical reactions in the future.

Exploring the Optimal Alloy for Nitrogen Activation by Combining Bayesian Optimization with Density Functional Theory Calculations.

Binary alloy catalysts have the potential to exhibit higher activity than monometallic catalysts in nitrogen activation reactions. However, owing to the multiple possible combinations of metal elements constituting binary alloys, an exhaustive search for the optimal combination is difficult. In this study, we searched for the optimal binary alloy catalyst for nitrogen activation reactions using a combination of Bayesian optimization and density functional theory calculations. The optimal alloy catalyst proposed by Bayesian optimization had a surface energy of ∼0.2 eV/Å2 and resulted in a low reaction heat for the dissociation of the N≡N bond. We demonstrated that the search for such binary alloy catalysts using Bayesian optimization is more efficient than random search.

Graph-theoretical exploration of the relation between conductivity and connectivity in heteroatom-containing single-molecule junctions.

In this study, we employ the Sachs graph theory to formulate the conduction properties of a single-molecular junction consisting of a molecule in which one carbon atom of an alternant hydrocarbon is replaced with a heteroatom. The derived formula includes odd and even powers of the adjacency matrix, unlike the graph of the parental structure. These powers correspond to odd- and even-length walks. Furthermore, because the heteroatom is represented as a self-loop of unit length in the graph, an odd number of passes of the self-loop will change the parity of the length of the walk. To confirm the aforementioned effects of heteroatoms on conduction in an actual sample, the conduction behavior of meta-connected molecular junctions consisting of a heterocyclic six-membered ring, whose conductive properties have already been experimentally determined, was analyzed based on the enumerated number of walks.

Photo-Induced Ring-Opening Reaction of Flav-3-en-2-ol Monitored by Time-Resolved Infrared Spectroscopy.

Photo-induced ring-opening reaction from flav-3-en-2-ol to 2-hydroxychalcone has been studied by time-resolved infrared (TR-IR) spectroscopy and quantum chemical calculations. A vibrational band due to the C═O stretching modes for intermediate species, enol forms of 2-hydroxychalcone in the electronic ground state, was observed at 1632 cm-1 in the TR-IR spectra after photoexcitation of flav-3-en-2-ol. We also found that the C═O stretching modes of the keto forms of 2-hydroxychalcone at 1664 cm-1 appeared immediately after photoexcitation and increased in intensity in synchronization with the depletion of the 1632 cm-1 band. Because the decay of the 1632 cm-1 band and the rise of the 1664 cm-1 band were fitted with bi-exponential model functions with common rate constants 0.5 and 11 µs-1, we propose that two kinds of enol form, single bond cis- (s-cis-) and trans- (s-trans-) enols, transformed into keto forms, cis-2-hydroxychalcone (Cc) and trans-2-hydroxychalcone (Ct), respectively. Quantum chemically calculated IR spectra of related species are consistent with the proposal. The observed temporal behavior of the TR-IR spectra indicates that there were reaction paths to the photogeneration of Cc and Ct within the time resolution of the TR-IR spectrometer (∼0.1 µs) in addition to the reaction paths via the enol forms of 2-hydroxychalcone.